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TROPICAL  CYCLONE  WIND  PROBABILITY  FORECASTING 
FOR  THE  NORTH  INDIAN  OCEAN  (WINDPIO) 


1 . 0  Introduction 

The  concepts  and  development  of  tropical  cyclone  wind 
probability  forecasting  for  the  North  Indian  Ocean  follow 
closely  those  originally  presented  by  Jarrell1  (1981)  for 
the  western  North  Pacific.  This  is  the  fourth  in  a  series 
of  wind  probability  forecasting  program  reports  covering 
different  ocean  basin  areas.  The  initial  report  on  the 
western  North  Pacific  was  followed  by  similar  reports  on  the 
western  North  Atlantic  (Jarrell2,  1981)  and  eastern  North 
Pacific  (Jarrell3,  1982).  This  report  and  the  two  pre¬ 
vious  are  based  upon  an  extension  of  the  concepts  developed 
for  tropical  cyclone  strike  probability  forecasting 
(Jarrell^,  1978).  As  these  concepts  have  been  previously 
presented  and  the  programs  are  presently  operational,  a 
redevelopment  of  the  concepts  will  not  be  presented  in  this 
report.  A  brief  summary  will  be  presented  as  a  background 


1Jarrell,  J.D.,  1981:  Tropical  Cyclone  Wind  Probability 
Forecasting  (WINDP),  NAVENVPREDRSCHFAC  Contractor  Report  CR 
81-03. 

2 

Jarrell,  J.D.,  1981:  Atlantic  Hurricane  Wind  Probability 
Forecasting  (WINDPA) ,  NAVENVPREDRSCHFAC  Contractor  Report  CR 
81-04. 

3Jarrell,  J.D.,  1982:  Tropical  Cyclone  Wind  Probability 
Forecasting  for  the  Eastern  North  Pacific  (EPWINDP), 
NAVENVPREDRSCHFAC  Contractor  Report  CR  82-06. 

^Jarrell,  J.D. ,  1978:  Tropical  Cyclone  Strike  Probability 
Forecasting,  NAVENVPREDRSCHFAC  Contractor  Report  CR  78-01. 


1 


for  the  new  reader.  Differences  between  the  development  of 
North  Indian  Ocean  tropical  cyclone  wind  probabilities  and 
that  of  its  closest  previous  counterpart,  the  model  for  the 
Atlantic,  is  the  crux  of  this  report  and  these  will  be 
described  in  detail. 

The  tropical  cyclone  wind  probability  forecasting 
model  generates  estimates  of  the  probability  of  30  and  50 
knot  winds  occurring  at  a  point  given  that  the  tropical 
cyclone  occupies  a  particular  position.  The  model  includes 
many  features  of  the  strike  probability  model,  which  is 
based  on  an  analysis  of  position  forecast  errors  to 
determine  the  probability  of  a  tropical  cyclone  occupying  a 
particular  geographic  position.  Wind  profile  errors  in  the 
forecast  of  maximum  wind  and  in  the  forecast  radius  of  30 
and  50  knot  winds  are  similarly  analyzed  to  determine  the  30 
and  50  knot  probabilities  of  occurrence. 

Jarrell^  (1978),  in  the  development  of  the  strike 
probability  model,  based  the  theory  of  strike  probability  on 
three  assumptions: 

1)  All  tropical  cyclone  forecasts  are  subject  to 
error; 

2)  Difficulty  of  forecast  and  size  of  forecast  error 
are  statistically  related;  and 

3 )  The  occurrence  of  errors  is  random  and  approx¬ 
imates  a  multimodal  bivariate  normal  probability 
distribution . 

Studies  by  independent  investigators  in  three  ocean 
basins  frequented  by  tropical  cyclones  verified  these 
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assumptions.  Nicklin5  (1977),  Thompson  et  al6  (1981), 
and  Crutcher  et  al7  (1982)  using  data  from  the  western 
Pacific,  eastern  Pacific  and  the  Atlantic,  respectively, 
also  developed  individual  methods  to  group  forecasts  into 
three  classes  of  forecast  difficulty.  A  general  relative 
classification  evolved  in  each  study  to  yield  forecast  error 
groups  of  below  average,  average,  and  above  average  errors. 
These  groups  were  also  referred  to  as  Class  I  (easy  fore¬ 
casts),  Class  II  (average  forecasts)  and  Class  III  (dif¬ 
ficult  forecasts)  in  previous  reports.  Each  investigator 
also  utilized  sufficient  statistical  data  and  prescribed 
sufficient  parameters  to  describe  the  bivariate  normal 
distributions  for  each  of  the  three  classes  and  for  fore¬ 
casts  of  24,  48,  and  72  hours. 

Following  the  work  of  Crutcher  et  al7  the  forecast 
errors  were  separated  into  groups  using  the  NORMIX  cluster¬ 
ing  model.  Unlike  the  other  basins,  only  two  significantly 
different  clusters  were  identified.  These  two  clusters  are 
shown  by  50%  probability  ellipses  in  Figure  1. 


D.S. ,  1977:  A  Statistical  Analysis  of  Western 

Pacific  Tropical  Cyclone  Forecast  Errors,  Naval  Postgraduate 
School,  M. S.  Thesis,  June. 

^Thompson ,  W.J.,  R. L.  Elsberry  and  R. G.  Read,  1981:  An 
Analysis  of  Eastern  North  Pacific  Tropical  Cyclone  Forecast 
Errors,  Monthly  Weather  Review,  V.  109,  pp.  1930-1938. 

7 

Crutcher,  H.L. ,  C.J.  Neumann  and  J.M.  Pelissier,  1982: 
Tropical  Storm  Forecast  Errors  and  the  Multimodal  Bivariate 
Normal  Distribution,  J.  Appl.  Meteor.,  V.  21,  pp.  978-987. 
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Figure  1.  Fifty  percent  error  probability 
ellipses  for  two  clusters  of  24-hr  forecast  errors 
identified  by  the  NORMIX  clustering  model  in  the 
North  Indian  Ocean.  Distances  are  in  nautical 
mi les  ( n  mi ) . 
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2.0 


Model  Description 


2.1  The  North  Indian  Ocean 

One  version  of  the  tropical  cyclone  wind  probability 
model  was  developed  for  the  Atlantic  Ocean  as  a  modification 
of  the  western  Pacific  wind  probability  program.  The  con¬ 
cepts  of  the  strike  portion  of  the  North  Indian  Ocean  model 
are  similar  to  those  developed  for  the  Atlantic.  The  North 
Indian  Ocean  model  uses  a  regression  equation  to  estimate 
the  probability  of  a  forecast  belonging  to  either  of  two 
forecast  classes  (associated  with  two  forecast  error  clus¬ 
ters)  and  the  Tsui,  et  al ,  model8  (used  in  the  Atlantic 
wind  probability  model)  to  determine  wind  profiles  and  to 
handle  asymmetrical  storm  distribution.  Thus  this  model 
uses  modified  versions  of  both  the  wind  probability  and 
strike  probability  concepts  used  in  the  Atlantic. 

The  Atlantic  wind  profiles  and  tropical  cyclone  asym¬ 
metry  are  derived  using  methods  produced  by  Tsui,  et  al ,  as 
indicated  above.  Tsui,  et  al,  using  wind  radius  data  from 
tropical  cyclone  warnings  over  a  12-year  period  (1966  to 
1977),  determined  that  the  profile  of  the  tangential  wind 
speed  along  the  radial  axis  was  exponential.  He  further  de¬ 
termined  that  maximum  wind  and  persistence  could  be  statis¬ 
tically  related  to  the  size  of  the  storm  and  that  the  asym¬ 
metric  shape  of  the  storm's  wind  pattern  could  be  correlated 


8Tsui ,  T.L. ,  L.R.  Brody  and  S.  Brand,  1982:  A  Technique 
for  Predicting  Surface  Wind  Distributions  of  Tropical 
Cyclones  in  the  Western  North  Pacific,  NAVENVPREDRSCHFAC 
Technical  Report  TR  82-05. 
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to  the  forward  speed  of  movement  of  the  storm.  A  simple 
empirical  relationship  was  derived  to  provide  an  estimation 
of  any  wind  radii: 


V/Vmax  =  exP  (-0.693R), 

where  V  is  the  wind  speed  of  interest,  is  the 

max 

maximum  wind  speed  of  the  storm,  and  R  is  a  ratio  of  the 
radius  associated  with  V  to  the  radius  associated  with 
one-half  of  Vmax  ,  respectively. 

Radii  are  measured  outward  from  the  radius  of  maximum  winds. 
Asymmetrical  storm  configuration  is  accommodated  in  the 
profile  by  an  empirical  adjustment  of  the  one-half  radius, 
rhalf'  dependent  on  bearing  relative  to  direction  of 
motion  and  translation  speed. 
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3 . 0  Testing  the  North  Indian  Ocean  Wind  Probability 

Program  (WINDPIO) 

The  methodology  used  in  comparing  WINDPIO  predicted 
values  against  observed  values  is  identical  to  that  used  in 
the  Atlantic.  An  array  of  30  points  in  the  North  Indian 
Ocean  was  selected  (Figure  2).  WINDPIO  values  for  the  30 
and  50  knots  were  calculated  at  12  hour  intervals  from  the 
effective  synoptic  time  of  the  Joint  Typhoon  Warning  Center 
( JTWC )  forecasts  for  the  three  1981  season  cyclones.  Since 
most  of  these  30  points  are  not  observing  stations,  actual 
verifying  winds  were  not  generally  available.  Consequently 
a  verifying  "warning  time"  probability  greater  than  50% 
constituted  a  verifying  strike. 

Tables  1,  2,  3  and  4  compare  the  expected  to  the 
observed  occurrences  of  30  and  50  knot  winds.  Predictions 
are  associated  with  percentage  groups  of  increasing  width, 
<1/2%,  1/2  to  1  1/2%,  ...,  etc.  Time  integrated  probabil¬ 
ities  were  verified  only  if  a  continuous  record  was  avail¬ 
able  over  the  entire  time  period.  Significance  of  the 
differences  between  the  expected  and  the  observed,  as  dis¬ 
cussed  in  previous  reports  (Jarrell9,  1981)  is  difficult 
to  assess,  but  using  a  "t"  test,  agreement  is  excellent. 

The  difference  between  the  expected  and  observed  occurrences 
were  never  at  a  significant  level  and,  although  the  sample 
is  extremely  small,  the  overall  results  are  considered 
statistically  sound. 


9Jarrell,  J.D.,  1981:  Atlantic  Strike  Probability 
Program,  NAVENVPREDRSCHFAC  Contractor  Report  CR  81-04. 
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I.u 


Figure  2.  Selected  test  points  C  30 )  in  the  North  Indian 
Ocean.  Included  are  seven  named  cities  of  interest  to  the 
U. S.  Navy. 


Table  1.  Instantaneous  Probabilities. 

50  kt  winds  -  Expected  versus  Observed. 


50  KT 

24  Hr  48  Hr  72  Hr 


A£P<B 

K 

E 
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K 
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K 

E 

< 
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312 
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41 
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2 

0 

14- 
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0 
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0 

0 

0 

ALL 
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1 

1 

330 

0 

0 

90 

0 

Table  2.  Instantaneous  Probabilities. 

30  kt  winds  -  Expected  versus  Observed. 

30  KT 

24  Hr  48  Hr  72  Hr 


A<P<B 

N 

E 

0 

K 

E 

0 

N 

E 

<  4% 

564 

0 

0 

260 
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0 

74 

0 

4-14 

49 

1 

1 

35 

0 

0 

11 

0 
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51 

1 

1 

35 

1 

1 

5 
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26 

1 

1 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 
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D 

0 

0 

0 

0 

0 

0 

0 

0 
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Table  3. 


Time  Integrated  Probabilities. 

50  kt  winds  -  Expected  versus  Observed. 

50  KT 

24  Hr  48  Hr  72  Hr 


A<P<B 

N 

E 

0 

K 

E 

0 

N 

E 

0 

610 

0 

0 

262 

0 

0 

67 

0 

0 

4-  14 

43 

0 

0 

43 

0 

0 

12 

0 

0 

14-  34 

32 

1 

1 

22 

1 

0 

11 

0 

0 

34-  74 

5 

0 

0 

3 

0 

1 

0 

0 

0 

74  -  154 

0 

0 

0 

0 

0 

0 

0 

0 

0 

134  -  314 

0 

0 

0 

0 

0 

0 

0 

0 

0 

314  -  634 

0 

0 

0 

0 

0 

0 

0 

0 

0 

>  634 

0 

0 

0 

0 

0 

0 

0 

0 

0 

ALL 

690 

1 

1 

330 

1 

1 

90 

0 

0 

Table  4.  Time  Integrated  Probabilities. 

30  kt  winds  -  Expected  versus  Observed. 

30  KT 

24  Hr  48  Hr  72  Hr 


A<P<B 

K 

E 

0 

N 

E 

0 

N 

E 

0 

<  4% 

493 

0 

0 

212 

0 

0 

55 

0 

0 

4  -  14 

80 

0 

1 

39 

0 

0 

8 

1 

0 

14-  34 

67 

1 

2 

56 

1 

1 

20 

0 

0 

34-  74 

37 

2 

2 

16 

1 

1 

7 

0 

1 

74  - 154 

13 

1 

1 

7 

1 

1 

0 

0 

0 

154  -  314 

0 

0 

0 

0 

0 

0 

0 

0 

0 

314  -  634 

0 

0 

0 

0 

0 

0 

0 

0 

0 

>  634 

0 

0 

0 

0 

0 

0 

0 

0 

0 

ALL 

690 

4 

6 

330 

3 

3 

90 

1 

1 
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4.0 


Operational  Products 


The  North  Indian  Ocean  wind  probability  program  will 

★  ... 

be  available  for  seven  preselected  points.  probabilities 
will  be  given  in  two  modes,  instantaneous  and  time  inte¬ 
grated,  and  at  0,  12,  24,  36,  48,  60  and  72  hours  after  the 
warning  time.  The  instantaneous  probability  will  be  the 
probability  at  the  stated  time  (i.e.,  12  hr)  and  the  time 
integrated  probability  will  be  summed  for  the  0  to  X  hour 
time  interval  for  an  estimate  of  the  probability  that  the 
event  will  be  observed  within  that  period  of  time. 

The  greatest  source  of  probable  error  for  the  WINDPIO 
program  will  be  erroneous  input  data.  An  internal  check  for 
unusual  motion  (expected  to  occur  only  5%  of  the  time  in 
nature)  will  be  made  and  suspect  motion  flagged.  The  user 
should  then  recheck  input  data  for  accuracy. 

When  the  forecast  track  approaches  a  land  mass,  the 
forecaster  should  be  aware  of  program  bias.  This  should  be 
minor  for  seaward  approach  to  low  coastal  areas  or  over 
smaller  islands.  However,  in  other  cases  land  influences 
will  appear  as  rapid  decreases  in  the  instantaneous  wind 
probabilities  (especially  50  kt  winds)  near  forecast  land¬ 
fall  time.  This  will  bias  probabilities  -  overstate  them 
for  inland  sites  and  understate  them  for  coastal  sites. 


*Locations:  Muscat,  Masirah  Island,  Karachi  and  Bombay  in 

the  Arabian  Sea,  and  Colombo,  Vishakhapatnam  and  Chittagong 
in  the  Bay  of  Bengal  (also  shown  in  Figure  1). 
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Time  integrated  probabilities  will  be  less  biased.  This 
problem  is  caused  by  wind  forecasts  being  influenced  by 
track  forecasts  where  landfall  is  concerned.  A  bad  track 
forecast  may  cause  a  bad  wind  forecast.  This  was  not 
accounted  for  either  in  development  nor  testing;  hence  the 
test  results  simulate  expected  actual  operational  results 
and  some  of  the  minor  disparities  between  expected  and 
observed  occurrences  no  doubt  stem  from  this. 

Figure  3  depicts  a  standard  output  of  WINDPIO  for  Cyclone 
31-81,  9  December  1981.  Probabilities  are  provided  at  12  hour 
intervals  for  strike,  50  and  30  kt  winds  in  both  instantaneous 
and  time  integrated  modes.  The  format  HHPIPS  where  HH  is  hours 
after  forecast  effective  time,  PI  and  PS  are  the  instantaneous 
and  time  integrated  probability  rounded  to  the  nearest  whole 
percent.  The  letters  "IN"  mean  insignificant  or  less  than  one 
half  percent  while  "THREAT  NIL"  means  all  probabilities  were 
less  than  one  percent. 
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STRIKE  AND  WIND  PROBABILITY  FORECAST 
CY  31-81  090800Z 


MUSCAT 

THREAT 

NIL 

MASIRAH  IS 

THREAT 

NIL 

KARACHI 

THREAT 

NIL 

BOMBAY 

50  KNOT 
30  KNOT 

00  ININ 
00ININ 
00  ININ 

121NIN 
12ININ 
12 ININ 

24  ININ 
24ININ 
24ININ 

36  ININ 
36  ININ 
36  ININ 

48ININ 

48ININ 

48IN01 

60  ININ 
60 ININ 
60  IN 02 

72ININ 
7  2 ININ 
72IN03 

COLOMBO 

50  KNOT 
30  KNOT 

00ININ 
00  ININ 
0  0 ININ 

12ININ 
12 ININ 
12ININ 

24 ININ 
24 ININ 
24 ININ 

36ININ 
36  ININ 
36ININ 

48ININ 
48ININ 
4 8 IN 01 

60  ININ 
60  ININ 
60IN02 

72ININ 

72ININ 

72IN03 

VISKAHPTNM 
50  KNOT 
30  KNOT 

00ININ 

000202 

002626 

120102 

120712 

123236 

240307 

240716 

242937 

360210 

360518 

362037 

480111 

480319 

481437 

600112 

600219 

601037 

720112 

720220 

720838 

CHITTAGONG 
50  KNOT 
30  KNOT 

00ININ 
00  ININ 
00ININ 

121NIN 

12ININ 

120204 

24IN01 

240103 

241014 

360204 

360409 

361921 

480308 

480513 

482225 

600211 

600516 

602127 

720213 

720517 

721928 

FOR  METEOROLOGISTS;  FORECAST  CONFIDENCE  TABLE 


TIME 

PROB 

DIST 

PROB 

DIST 

PROB 

12  HR 

15 

50 

15 

75 

70 

2  4HR 

24 

100 

22 

150 

54 

48HR 

33 

200 

26 

300 

41 

72HR 

38 

300 

27 

450 

35 

PROBABILITIES  BASED  ON  FOLLOWING  FORECAST 

001640866075  121720869070  241800872065  481980892045  722080922030 


Figure  3.  Sample  WINDPIO  output  for  cyclone  30-81,  9  December 
1981.  See  text  for  format.  This  product  is  based  on  an 
official  forecast  issued  by  the  Joint  Typhoon  Warning 
Center,  Guam. 
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5 . 0  Summary 


The  wind  probability  model  for  the  North  Indian  Ocean 
is  largely  based  on  the  wind  probability  program  for  the 
Atlantic.  Many  features  of  wind  probability  programs  for 
the  eastern  and  western  North  Pacific  have  also  been 
incorporated.  All  of  these  programs  are  operational  or  soon 
to  be  operational.  Test  results  of  the  North  Indian  Ocean 
wind  probability  program  (WINDPIO)  demonstrated  excellent 
agreement  between  expected  and  observed  results.  Operatioal 
products  are  available  both  in  message  form  as  described 
herein  and  in  the  Navy  Environmental  Display  System  (NEDS) 
as  plotted  charts. 
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